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Abstract Accretion onto supermassive black holes produces both the dramatic 
phenomena associated with active galactic nuclei and the underwhelm- 
ing displays seen in the Galactic Center and most other nearby galax- 
ies. I review selected aspects of the current theoretical understanding 
of black hole accretion, emphasizing the role of magnetohydrodynamic 
turbulence and gravitational instabilities in driving the actual accre- 
tion and the importance of the efficacy of cooling in determining the 
structure and observational appearance of the accretion flow. Ongo- 
ing investigations into the dynamics of the plunging region, the origin 
of variability in the accretion process, and the evolution of warped, 
twisted, or eccentric disks are summarized. 



3.1 Introduction 

Most galactic nuclei are now believed to harbor supermassive black 
holes. These black holes are all accreting gas — at a minimum from the 
interstellar medium proximate to the event horizon, and in some cases 
from dense disks of gas in orbit around the hole. The observational 
signatures of this accretion, however, differ dramatically from galaxy to 
galaxy. The black hole in the Galactic Center, for example, has an X- 
ray luminosity in quiescence of only Lx — 2 x 10'^^ ergs s^^, reaching 
Lx ~ 10^^ ergs s~^ during flaring states (Baganoff et al. 2001, 2003; 
Goldwurm et al. 2003). Many nearby elliptical galaxies also show nu- 
clear emission that is much weaker than might be expected on the basis 
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of simple estimates of the black hole accretion rate (Fabian & Rees 1995; 
Di Matteo et al. 2000; Loewenstein et al. 2001). In stark contrast to 
these feeble displays, accretion onto black holes in quasars powers the 
most luminous steady sources in the universe. At a minimum, a the- 
ory of accretion needs to account for this dichotomy, and to explain at 
least some of the many secondary phenomena (jets, outflows, variability, 
etc.) associated with Active Galactic Nuclei (AGN). More ambitiously, 
one might hope to understand the role of accretion in actually forming 
super massive black holes. Black hole formation is a difficult theoretical 
problem which, although currently untroubled by direct observations, is 
receiving increasing attention. 

In this chapter, I discuss selected aspects of the theory of accretion 
onto supermassive black holes, with an emphasis on the physical pro- 
cesses that drive accretion and determine the qualitative properties of 
the flow. In §3.2, I outline the mechanisms which can lead to angular 
momentum transport within an accretion flow, thereby allowing rota- 
tionally supported gas to flow inward and liberate gravitational potential 
energy. Close to the black hole, turbulence driven by magnetohydrody- 
namic instabilities is probably the dominant mechanism for transport. 
Further out, other processes, such as gravitational instabilities, are likely 
to assume that role. Angular momentum transport, however, is only part 
of the story. Although central to our understanding of accretion, knowl- 
edge of its origin does not suffice to explain why the Galactic Center 
looks nothing like a powerful AGN. For that we need to consider the 
distinction between geometrically thin accretion disks, in which the gas 
can radiate efficiently and cool to sub-virial temperatures, and hot thick 
disks, which are radiatively inefficient and vulnerable to rapid mass loss. 
These questions are addressed in §3.3. Subsequent sections examine the 
status of several open questions in the study of black hole accretion, 
including the dynamics of gas executing its final plunge into the hole, 
the origin of variability, and the evolution of disks that are warped or 
eccentric. 

3.2 Angular Momentum Transport 

In almost all circumstances, gas in the nuclear regions of galaxies 
has far too much angular momentum to be swallowed directly by the 
black hole. Understanding the mechanisms that can lead to angular 
momentum transport in rotating flows is therefore the central problem 
in the theory of accretion onto supermassive black holes. At scales of 
10 — 100 pc or larger, this is a problem in galactic dynamics. Following 
galactic mergers, for example, gravitational torques from bars and other 



transient structures can efficiently funnel gas into the nuclear regions 
(Shlosman, Prank, & Begelman 1989; Barnes & Hernquist 1991; Hern- 
quist & Mihos 1995). Once within the black hole's sphere of influence, 
however, which for a galaxy of central velocity dispersion a and black 
hole mass M, extends out to 

GM, 

r, = (3.1) 

these galactic mechanisms become less efficient. Within the sphere of 
influence, which reaches 10 pc for a black hole of mass M, = 10^ M©, 
the black hole's gravity overwhelms that of the host galaxy and domi- 
nates the potential, unless a comparably large mass of gas has managed 
to accumulate at such small radii. Timescale arguments suggest that 
gas at r, probably cannot trickle down through an accretion disk all 
the way to the vicinity of the black hole (e.g., Shlosman, Begelman, & 
Frank 1990). Disks probably do form, however, at smaller radii of order 
0.1 pc, where maser emission in the nucleus of NCG 4258 has an un- 
mistakably disk- like geometry (Miyoshi ct al. 1995). At such radii, the 
specific angular momentum of a Keplerian disk exceeds that needed for 
direct capture by the black hole by a factor of order 10^. If that gas is 
to accrete, 99% of the angular momentum must either be redistributed 
within the disk — angular momentum transport — or lost entirely from the 
system. We will consider the known mechanisms which can accomplish 
this feat. 

In the inner disk, turbulence driven by mxignetic instabilities is a po- 
tent source of angular momentum transport. Although magnetohydro- 
dynamic (MHD) disk turbulence has long been implicated in disk an- 
gular momentum transport and is mentioned prominently by Shakura 
& Sunyacv (1973), current confidence in this conclusion rests on two 
more recent developments. First, Balbus & Hawley (1991) demonstrated 
that the introduction of a weak magnetic field renders accretion flows 
linearly unstable to a powerful local instability. Subsequent numerical 
simulations (Hawley, Gammie, & Balbus 1995; Brandenburg et al. 1995; 
Matsumoto &: Tajima 1995; Stone et al. 1996) showed that the instabil- 
ity rapidly develops into sustained turbulence, which transports angular 
momentum outward at a rate that is consistent with observational con- 
straints derived from studies of accretion in mass transfer binaries (e.g., 
Pringle, Verbunt, & Wade 1986; Cannizzo 1993; Osaki 1996; Hameury et 
al. 1998). For a recent review of the role of MHD turbulence in accretion 
disks, see Balbus (2003). 

In the absence of magnetic fields, a differentially rotating disk with 
angular velocity ^{r) is linearly stable to axisymmetric perturbations 



according to the Rayleigh criterion if 




(3.2) 



i.e., if the specific angular momentum /(r) of the flow is an increasing 
function of radius. For a geometrically thin disk orbiting a point mass 
M, the angular velocity is Keplerian, 



the specific angular momentum l{r) oc ^/r, and the disk is hydrodynam- 
ically stable. Numerical simulations support this conclusion (Balbus, 
Hawley, & Stone 1996). Although real disks — especially geometrically 
thick flows with significant pressure support — can have angular velocity 
profiles that differ from the simple Keplerian form, they too are invari- 
ably stable by the Rayleigh criterion. 

Matters are drastically different if the disk contains a magnetic field. 
Analytic studies have shown that a weak magnetic field destabilizes as- 
trophysical disks, provided that 



a condition which is almost always satisfied in real disks. This mag- 
netorotational instability (MRI) exists regardless of the initial magnetic 
field configuration (Velikhov 1959; Chandrasekhar 1960; Balbus & Haw- 
ley 1991, 1992; Terquem & Papaloizou 1995; Ogilvie & Pringle 1996; 
Curry & Pudritz 1996), though the growth rates vary depending upon 
the magnetic field geometry, being fastest for vertical fields. 

Proving the existence of the MRI in the general case requires a mod- 
erately involved calculation, which can be found in the comprehensive 
review by Balbus &: Hawley (1998). Analyses of much simpler systems, 
however, reveal most of the important physics. Here I follow closely the 
treatment of Balbus &: Hawley (2000) . Consider a fluid element orbiting 
in a disk with a central gravitational potential $(r). We will assume 
that pressure forces are negligible, as they will be, provided that per- 
turbations to the disk velocity field are highly subsonic. In cylindrical 
polar coordinates (r, ^), the equations of motion then read 




(3.3) 



r — r<j)'^ = 
rcf) + 2f(f) = 



(3.5) 



where the dots denote derivatives with respect to time, and fr and are 
forces that wc will specify shortly. We now concentrate our attention 
on a small patch of the disk at radius tq that is corotating with the 
overall orbital motion at angular velocity We define a local Cartesian 
coordinate system {x, y) via 



(3.6) 
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and substitute these expressions into equation (3.5) above. Discarding 
quadratic terms, the result is 



X — 2Qy 
y + 2Qx 
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dlnr 



(3.8) 



These equations describe the epicyclic motion of pressureless fluid per- 
turbed from an initially circular orbit. 

If the disk contains a weak vertical magnetic field, perturbations to 
the fluid in the plane of the disk will be opposed by magnetic tension 
forces generated by the bending of the field lines. Considering in-plane 
perturbations varying with height z and time t as e^^'^^~^^) ^ the magnetic 
tension force is f = — (A;t;yi)^s, where s is the displacement vector and 
VA = ^/i?|/47^p is the Alfven speed. Using this expression for fx and /y, 
and assuming the e*'^* time dependence, equation (3.8) becomes 



-oJ^x — 2iu^y 
-oj^y + 2iuQ,x 



dlnr 

-{kvAfy- 



{kvA^x , 



(3.9) 
(3.10) 



Combining these equations yields a dispersion relation that is a quadratic 
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= 0. 
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The system is unstable if < 0, which occurs when the third term in 
the dispersion relation is negative. Instability therefore occurs if 



< 0. 



(3.12) 



For a sufficiently weak field {va 0), or for long enough wavelength 
perturbations (A; — 0), we then obtain the aforementioned criterion for 
a disk to be unstable to the MRI, namely 



dJ^2 



< 0. 



(3.13) 



dlnr 



In a real disk, the longest wavelength perturbation in the vertical direc- 
tion will be of order the scale height h of the disk. There will therefore be 
instability, provided that the magnetic field is weaker than some thresh- 
old value -Bmax- 

The physical origin for this instability is fairly straightforward. Adopt- 
ing the same configuration that we have just analyzed, consider the effect 
of perturbing a weak vertical magnetic field threading an otherwise uni- 
form disk. If the field remains frozen into the plasma, then field lines 
which connect adjacent annuli in the disk will be sheared by the differ- 
ential rotation of the disk into a trailing spiral pattern. Provided that 
the field is weak enough, magnetic tension is inadequate to snap the 
field lines back to the vertical. What tension there is, however, acts to 
reduce the angular momentum of the inner fluid element and boost that 
of the outer fluid element, providing angular momentum transport in 
the outward sense that is required to drive accretion. 

The pervasive nature of the MRI, taken together with its rapid linear 
growth rate — <^max can be as large as 317/4, meaning that growth of 
magnetic field energy occurs on a dynamical timescale — implies that 
well-ionized astrophysical disks will be unstable. Analytic studies are 
unable, however, to investigate the most interesting consequence of the 
instabilities, namely, what happens when the MRI reaches a saturated 
state and turbulence has developed throughout the flow? Quantities we 
would like to determine in this state include the rep component of the 
stress tensor. 



which includes both fluid (Reynolds) and magnetic (Maxwell) stresses. 
In this expression, dv^ = — r^l is the fluctuation in the azimuthal 
velocity, v^j- and va^ are Alfven velocities in the radial and azimuthal 
directions, respectively, and the angle brackets denote a density weighted 
average over height. If Wr(j) can be measured, in practice from a numer- 
ical MHD simulation, then it can be used to estimate the a parameter 
introduced by Shakura & Sunyaev (1973). They suggested the scaling 
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where Cg is the sound speed and a is a dimensionless constant or pa- 
rameter that measures the efficiency of angular momentum transport in 
disks. 

Numerical simulations of the non-linear development of the MRI have 
now been performed in both local and global geometries. Local simula- 
tions (Hawley, Gammie, & Balbus 1995, 1996; Brandenburg et al. 1995; 
Matsumoto & Tajima 1995; Stone et al. 1996; Miller & Stone 2000) 
follow the evolution of the MRI in a small patch of disk, with periodic 
boundary conditions in <p and (in a modified form) r. Global simula- 
tions (Armitagc 1998; Hawley 2000) have the advantage of being able 
to study large scale magnetic fields and are essential for the investiga- 
tion of geometrically thick {h ~ r) flows (e.g., Hawley, Balbus, & Stone 
2001). All such simulations have known limitations. The physical sepa- 
ration between the largest relevant scale (that of the whole disk, or the 
wavelength of the most important MRI modes) and the dissipative scale 
(either viscous or resistive) is almost always too large to resolve numer- 
ically. Moreover, many calculations rely on purely numerical effects at 
the grid scale to provide dissipation. Discussions of the effect of these 
limitations can be found, for example, in Balbus &; Hawley (1998) or 
Schekochihin et al. (2004). Nevertheless, it is encouraging that a vari- 
ety of simulations, using different numerical techniques, agree in several 
critical aspects. The MRI rapidly leads — on a timescale of just a few 
orbital periods — to a state of MHD disk turbulence in which angular mo- 
mentum is transported outward. This is sustained, despite the presence 
of dissipation (often at unphysically large scales) within the numerical 
codes. The magnetic field energy remains smaller than the thermal en- 
ergy in the disk. Finally, the Maxwell stress dominates by a large margin 
over the Reynolds stress, though the latter also provides some outward 
transport. Values of a between 10~^ and 10~^ are representative of the 
bottom line of most simulations. 

The maintenance of a disk magnetic field, despite the presence of 
dissipation, implies that turbulence driven by the MRI sustains a disk 
dynamo, just as fiuid motions in stars and planets generate their own 
self-sustaining magnetic fields (Parker 1955; Proctor &; Gilbert 1994; 
Glatzmaicr & Roberts 1995). The accretion disk case, however, is inter- 
estingly different from these better-observed systems. In a disk (at least 
according to current belief) there would be no turbulent motions in the 
absence of a magnetic field, which rather must generate the turbulent 
velocity field needed for its own amplification. This bootstrapping qual- 
ity makes it harder to apply knowledge gleaned from the long history of 
work on solar and planetary dynamos to the disk case. 



How do the values of a obtained from numerical simulations compare 
with those inferred from observations? The best constraints come from 
modeling the episodic accretion events seen in the dwarf novae subclass 
of cataclysmic variables. In the widely accepted thermal disk instability 
model for dwarf nova outbursts (Meyer & Meyer-Hofmeister 1981; Mi- 
neshige & Osaki 1983; Faulkner, Lin, & Papaloizou 1983), the eruptive 
behavior results from the inability of the disk to reach a steady-state at 
radii where hydrogen is partially ionized. Under these conditions, there 
is no thermally stable vertical structure for a range of accretion rates, 
resulting in a limit cycle in which the disk alternates between a cool 
mass-accumulating phase and a hot draining phase. The timescales of 
the outbursts are directly related to the values of a in the outburst and 
quiescent states. Different authors (e.g., Pringle et al. 1986; Cannizzo 
1993; Osaki 1996; Hamcury et al. 1998) derive typical values of a in the 
outburst state of ~ 0.1, with quiescent values a factor of a few lower. 
This is in reasonably encouraging agreement with the values predicted 
from MHD simulations. Moreover, as I will discuss later in the context 
of AGN disks, the difference in a between the outburst and quiescent 
states could be caused by a suppression of the MRI when the disk is 
cool and the magnetic field is imperfectly coupled to the gas (Gammie 
& Menou 1998). 

Of course, in many circumstances of interest, it is either inconve- 
nient or unfeasible to determine the structure and evolution of disks 
via three-dimensional MHD simulations. In particular, the evolution of 
AGN disks across the full range of radii spanned by real disks cannot 
be simulated directly. This raises the following questions: under what 
circumstances are vertically-integrated one-dimensional disk models, or 
two-dimensional simulations using a parameterized viscosity, good ap- 
proximations to the real disk physics? 

For geometrically thin disks [hjr <^ 1), the standard time-dependent 
one-dimensional treatment divides the problem into two parts (e.g., 
Frank, King, & Raine 2002). First, one solves for the vertical structure 
of an isolated annulus of the disk using methods and equations analogous 
to those of stellar structure. The vertical structure calculation yields the 
functional dependence of the kinematic viscosity v on the surface den- 
sity S, angular velocity O, and a. Armed with this knowledge, we can 
then solve a one-dimensional diffusion equation for the time-dependence 
of the disk, 



subject to appropriate boundary conditions at the inner and outer edges 
of the disk. If thermal fronts are present (or develop) within the disk. 




(3.16) 



this equation for S needs to be supplemented with an explicit equation 
for the time-dependence of the central temperature Tc (e.g., Cannizzo 
1993). 

Provided that we restrict our attention to the axisymmetric evolution 
of geometrically thin disks, one-dimensional disk models constructed in 
this manner probably provide a reasonable approximation to the evolu- 
tion of MRI-active disks. At the most basic level, the MRI is a local 
instability, and, as such, is in principle amenable to an a-type descrip- 
tion (Balbus & Papaloizou 1999). For the vertical structure part of the 
calculation, the main question is whether the distribution of the stress 
with height follows the simple scalings (stress proportional to pressure) 
normally employed as a generalization of the a ansatz (Eq. 3.15). In fact, 
numerical simulations of the MRI suggest that the stress may be roughly 
constant within the body of the disk (up to ~ 2h), with a significant 
fraction (perhaps a quarter) of the magnetic energy escaping buoyantly 
from the disk to form a hot, strongly magnetized corona (Miller & Stone 
2000). These simulation results differ from the simplest assumptions 
used in modeling vertical structure, but they do not invalidate the basic 
approach. Indeed, it ought to be possible to compute vertical structure 
models — and vital observational properties of the disk, such as the lo- 
cal spectrum (e.g., Hubeny & Hubeny 1997; Hubeny ct al. 2001) — using 
the computational results as a guide to where in the disk dissipation of 
energy occurs. This has not yet been done. 

The second fundamental assumption of one-dimensional models is the 
diffusive evolution of the surface density predicted by Equation 3.16. 
This, too, is broadly compatible with what is known about the MRI. Fig- 
ures 3.1 and 3.2 show the predicted evolution of a ring of gas in a Kcplc- 
rian potential, computed directly from a global three-dimensional MHD 
simulation using methods described by Armitage (1998). Although large 
scale non-axisymmetric structure is evident in Figure 3.1, the diffusive 
evolution of the surface density that is predicted by the one-dimensional 
theory is nonetheless reproduced. 

The preceding discussion has emphasized that one-dimensional disk 
models based on the a formalism can provide a reasonable description 
of the evolution of thin disks, especially if insights from numerical simu- 
lations of the MRI are used to guide the detailed modeling. Simulations, 
however, also reveal circumstances under which a models fail to describe 
disks driven by magnetically induced turbTilcncc. Papaloizou & Nelson 
(2003) and Winters, Balbus, & Hawley (2003a) have shown that av- 
eraging the magnetically induced stress over tens or even hundreds of 
orbital times is required to yield stable estimates of the average stress. 
In the context of a one-dimensional a disk model, this implies that a 



Figure 3.1. Evolution of a narrow ring of gas under the influence of angular mo- 
mentum transport, computed using a three-dimensional global MHD simulation (Ar- 
mitage, unpublished). Panels show the surface density E(r, 0) at different times t, 
where t is measured in units of the orbital period at the initial center of the ring. The 
simulation is isothermal, with Cs/v^ = 0.2 at r = ro, the initial center of the ring. The 
MRI was triggered by threading an initially stable and axisymmetric ring of gas with 
a weak vertical magnetic field. A cylindrical disk approximation was used (i.e., no 
vertical stratification), and the evolution was computed using the ZEUS code (Stone 
& Norman 1992) with 200 x 250 x 40 grid points in r, 0, and z, respectively. As in 
the one- dimensional viscous disk solution to this problem, given by Lynden-Bell & 
Pringle (1974), the ring spreads out, with the mass flowing inward, while a fraction of 
the mass at large radii absorbs the angular momentum. In addition, however, there is 
obvious non-axisymmetric structure in the form of trailing spiral waves and relatively 
long-lived fluctuations in the magnitude of the derived stress. 
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Figure 3.2. Evolution of the azimuthally-avcragcd surface density profile (in arbi- 
trary units) from an MHD simulation of the spreading of a narrow ring of gas. Curves 
correspond to the images shown in Fig. 3.1 {t = 0, 2, 4, and 6 orbits at the initial 
location of the ring). The simulated disk profile remains smooth. 

cannot be regarded as a constant over shorter timescales. More gener- 
ally, any attempt to extend a models to treat two-dimensional flows (for 
example, to study eccentric disks or disks strongly perturbed by binary 
companions) is dangerous and not to be recommended. Such generaliza- 
tions frequently rest on the assumption that all the components of the 
stress tensor in a turbulent MHD disk flow behave in the same way as 
a Navier-Stokcs viscosity. There is scant reason to believe that this is 
true, and indeed the limited numerical evidence collected to date sug- 
gests that the multidimensional evolution of magnetically driven disks 
differs significantly from that of their viscous counterparts (Nelson & 
Papaloizou 2003; Winters, Balbus, k Hawley 2003b). 

So far, we have assumed that the disk orbits in the point mass po- 
tential of the black hole, ignoring the self-gravity of the disk itself. This 
approximation is generally valid in the inner regions of the disk and will 
obviously fail if the disk mass becomes comparable to the mass of the 
black hole (perhaps during formation of the hole via accretion). Even if 
the disk mass is relatively modest, however, self-gravity can still become 




important at large enough radii. To demonstrate this, we make use of 
the a-prescription and write the kinematic viscosity v in the form 



For a steady disk, away from the boundaries (e.g., Pringle 1981), 

M 

i.S = -, (3.18) 

which, together with Equation 3.17, specifies S(r) in terms of tlie ac- 
cretion rate, sound speed, and a parameter. The stabihty of a nearly 
Keplerian gas disk against axisymmetric perturbations is measured by 
the Toomre (1964) Q parameter, 

with instabilities developing as Q 1 from above. Making use of Equa- 
tions 3.17 and 3.18, we find that for a steady disk, 

3 

Q OC ^ . (3.20) 

M 

Since the temperature and sound speed in the disk are decreasing func- 
tions of radius (at least if the heating is provided by accretion and/or 
irradiation from the central regions of the disk), the outer parts of a 
steady-state disk will inevitably reach the Q ~ 1 threshold for self- 
gravity to become important, provided that the disk is large enough. 
Indeed, models of AGN disks suggest that self-gravity typically sets in at 
r ~ 10~^ pc, which is only ~ 10^ Schwarzschild radii for a 10^ Mq black 
hole (Clarke 1988; Kumar 1999; Goodman 2003). Interestingly, this ra- 
dius is one to two orders of magnitude smaller than the radii from which 
maser emission has been observed in NGC4258 (Miyoshi et al. 1995) and 
NGC1068 (Grccnhill & Gwinn 1997). This raises the immediate ques- 
tion: is the masing disk observed in these galaxies a roughly homogenous 
structure, or is it rather composed of discrete clumps (Shlosman et al. 
1990; Kumar 1999; Kartje, Koiiigl, & Elitzur 1999)? 

Once Q « 1 has been attained, self-gravity triggers the development of 
spiral structure in the disk. Gravitational torques then act to transport 
angular momentum outward, in a manner superficially resembling the 
evolution of an ordinary viscous disk described by an a model (Laughlin 
&; Bodenheimer 1994). Such a correspondence rests on shaky theoretical 
foundations (Balbus &; Papaloizou 1999), because the non-local nature 




Figure 3.3. Surface density from global smooth particle hydrodynamics simulations 
of gravitationally unstable disks. The disk in this simulation had a mass equal to one 
tenth of that of the central object. In the simulation depicted in the left-hand panel, 
the cooling time tcooi = 5 Q,~^ . Under these conditions, a balance between heating 
and cooling is attainable, and the disk develops a stable spiral structure that acts 
to transport angular momentum outward. A shorter cooling time of tcooi ~ 3 
(right-hand panel) leads to prompt fragmentation of the disk into multiple bound 
objects. Note that the thickness of the simulated disk exceeded that appropriate for 
AGN. Figure adapted from Rice et al. (2003; their Figs. 2 and 4). 

of gravitational angular momentum transport immediately violates the 
local property of a models (Eq. 3.15). That said, examining the local 
properties of self-gravitating disks proves to be useful in delineating the 
circumstances under which a self-gravitating disk will fragment, rather 
than stably transport angular momentum (Lodato & Rice 2004). If, at 
a particular radius, the disk can cool via radiative losses on a timescale 
tcooi) then Gammie (2001) showed that the boundary between stable 
transport and fragmentation lies at tcooi = ^crit — 3 For tcooi > ^critj 
a stable state can be reached in which Q ~ 1 and angular momentum is 
transported with an effective a of 

-1 

(3.21) 

where 7 is a two-dimensional adiabatic index (Gammie 2001). As icooi 
decreases, turbulence in the disk (measured by the effective a) has to be- 
come increasingly violent in order to balance the radiative losses, and the 
overdensities of spiral structures increase. Eventually, for tcooi < ^crit, 
the disk fragments into bound objects. Both local (Gammie 2001; John- 
son & Gammie 2003) and global (Rice et al. 2003; Fig. 3.3) simulations 
broadly support these analytic results. 
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Application of these ideas to the outer regions of accretion disks in 
AGN is complicated by the dominant role that irradiation plays in set- 
ting the thermal conditions in the outer disk. In particular, fluctuations 
in the disk temperature (and, thus, Q) will occur as the illumination 
from the central source varies, with the disk being most vulnerable 
to fragmentation when the source is in its lowest state. Furthermore, 
there is no reason to expect that the disk at large radii is able to at- 
tain the steady-state described by Equation 3.18. The most plausible 
scenario is that the inner non-self-gravitating disk is surrounded by a 
region in which self-gravity transports angular momentum in the stable, 
non-fragmenting regime. Beyond that, at radii of order a pc or less, 
any disk formed from inflowing or outflowing gas is unstable to frag- 
mentation (Goodman 2003). If correct, this suggests that AGN need to 
be resupplied with gas of low specific angular momentum (compared to 
that of a circular orbit at the sphere of influence) in order to prevent 
fragmentation occurring in the outer disk. More speculatively, if frag- 
mentation does occur, the edge of the accretion disk in AGN could be a 
birth place for massive stars that would ultimately form stellar or inter- 
mediate mass black holes. If this process had occurred in the Galactic 
Center in the recent ('^lO Myr) past, it could leave a distinctive signa- 
ture in the form of disk-like kinematics for some subset of massive stars 
(Levin &; Beloborodov 2003). Moreover, if some fraction of the stars 
ended up embedded in the inner disk, they would migrate inward on 
at most a viscous timescale. Gravitational waves from the final inspiral 
of the stellar remnants with the supermassive black hole are potentially 
detectable with the Laser Interferometer Space Antenna {LISA; Levin 
2004; Goodman k Tan 2004). 

In broad outline then, the MRI is likely to provide the main mecha- 
nism of angular momentum transport in the inner disk, while self-gravity 
is potentially important at relatively large radii. A more subtle ques- 
tion is whether there is a region of the disk in which both processes are 
operative, or, conversely, a zone in which neither drives efficient angu- 
lar momentum transport. This depends upon whether the MRI is able 
to operate in the cool outer regions of the disk, where the gas has a 
low ionization fraction and the flux-freezing assumption of ideal MHD 
breaks down. For AGN, a quantitative discussion of the effect of the 
low temperature on angular momentum transport processes has been 
given by Menou Sz Quataert (2001). For a disk with sound speed Cg and 
resistivity r/, the magnetic Reynolds number is 



(3.22) 



If ReM is too low, then diffusion of magnetic fields due to the finite 
conductivity can suppress the MRI (e.g., Matsumoto &; Tajima 1995; 
Gammie 1996; Wardlc 1999; Sano & Stone 2002). The critical value of 
Rcm (which may also need to take account of the Hall effect) is probably 
of order unity or somewhat larger. Identical reasoning suggests that 
ambipolar diffusion can also shut off the MRI and prevent MHD disk 
turbulence from developing (Mcnou & Quataert 2001). 

Applying these ideas to AGN disks, Menou & Quataert (2001) found 
that the ionization state of the outer disk — and hence the likelihood of 
the MRI being operative — was highly uncertain, depending upon (among 
other things) whether the outer disk was exposed to irradiation from the 
central source. If the MRI does shut off at radii where Q » 1, then an 
intermediate zone in the disk could be both magnetically and gravita- 
tionally inactive. Such a "dead zone" would act to starve the inner disk 
of resupply from larger radii, until enough matter built up to initiate 
self-gravity or magnetorotational activity. Complicated and poorly un- 
derstood time-dependent accretion would probably result (e.g., Gammie 
1999b; Armitage, Livio, & Pringle 2001). Alternatively, if the regions of 
MHD turbulence and gravitational instability overlap, then there could 
be non-trivial interactions between the two angular momentum trans- 
port mechanisms. Planet formation simulations (Winters et al. 2003b) 
show that external gravitational torques can modify the properties of 
MHD disk turbulence, and the first simulations of self-gravitating mag- 
netized disks (Fromang et al. 2004) likewise show that magnetic and 
gravitational torques are not simply additive. In the calculations of Fro- 
mang et al. (2004), the presence of MHD turbulence in a self-gravitating 
disk significantly reduced the rate of gravitational transport of angular 
momentum. 

Finally, it remains possible that accretion is not always driven by 
angular momentum redistribution within the disk at all, but rather by 
angular momentum loss in an outflow. Outflows and jets are certainly 
present in some AGN, and if magnetically driven (Blandford &; Payne 
1982), they can potentially remove alloi the angular momentum needed 
to allow accretion, while ejecting only a modest fraction of the mass. 
Models in which outflows arc entirely responsible for driving accretion 
are not popular, in part because of skepticism that the required mag- 
netically launched outflows are present across a sufficiently broad range 
of disk radii. In particular, although there is near-universal agreement 
that magnetic fields are implicated in the formation of rclativistic jets, 
the slower outflows that originate further away from the black hole may 
well be partially or entirely accelerated by radiation pressure acting on 
resonance lines (Vitello & Shlosman 1988; Murray et al. 1995; Proga 



2003). Moreover, the stability of accretion arising from magnetic out- 
flow torques is rather doubtful (Lubow, Papaloizou, & Pringle 1994; Cao 
& Spruit 2002; but see also Konigl & Wardle 1996 and Campbell 2001). 
Nonetheless, the possibility that angular momentum loss via MHD out- 
flows dominates over internal stresses in restricted regions of AGN disks 
remains entirely viable. It constitutes a significant source of uncertainty 
in our understanding of the structure of the disk and when considering 
possible mechanisms for large-scale variability. 

3.3 The Effect of the Accretion Rate 

The accretion rate onto a black hole is most usefully measured as 
a fraction of the accretion rate at which radiation pressure becomes 
important for the dynamics of the flow. For a spherically symmetric 
flow, radiation pressure acting on free electrons balances gravity at the 
Eddington luminosity, 

LEaa = ^^M, (3.23) 

where gt is the Thomson scattering cross-section. If the luminosity 
is generated by accretion with radiative efficiency e, then L = tMc? 
and the above expression defines an Eddington accretion rate. Taking 
e = 0.1, 

^Edd = 2.2 f-^) Mo yr-i . (3.24) 
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The accretion rate can then be expressed in dimensionless form as a 
fraction m = M/MEdd of the Eddington rate. Current theoretical un- 
derstanding identifies rh as the most important parameter controlling the 
structure and observational appearance of black hole accretion flows. We 
will consider in turn the three cases: m <C 1, m ~ 1, and m 3> 1. These 
cases, by definition, differ in the importance of radiation pressure for the 
dynamics of the flow. Equally important, the ability of the flow to cool 
also varies systematically with rh. 

For m ^ 1, the critical question is whether the inflowing gas can 
radiate the gravitational potential energy liberated by accretion. Most 
theoretical work suggests that the vast majority of the energy dissipated 
by MHD turbulence in a disk goes initially into heating the ions rather 
than the electrons (Quataert 1998; Blackman 1999; Quataert Sz Gruzi- 
nov 1999; Medvedev 2000), at least as long as the magnetic field remains 
sub-equipartition. Sub-equipartition fields are the normal outcome of 
the non-linear phase of the MRI in MHD simulations of accretion fiows 
in the fluid limit (for a discussion of how the MRI itself is altered if 



the plasma is coUisionless, see Quataert, Borland, & Hammett 2002). 
Radiative losses from bremsstrahlung radiation, synchrotron radiation, 
and Comptonization of soft photons, on the other hand, primarily take 
energy from the electrons. In this situation, where the ions are heated 
simultaneously with the electrons being cooled, the extent of coupling 
between the two fluids obviously determines the state of the plasma. If 
the dominant coupling process is ordinary Coulomb collisions, then the 
timescale for thermal equilibrium to be established between the ions and 
the electrons is (Spitzer 1962; Ichimaru 1977) 



Here, Tg is the electron temperature, k is the Boltzmann constant, rup 
is the proton (ion) mass, n is the number density of ions, and rrie is 
the electron mass. The Coulomb logarithm has a value of In A ~ 20. 
When the accretion rate is low, n is small, and this, coupled with the 
high temperature near the black hole (of the ions generally, and the 
electrons, as well, if we initially postulate a onc-tcmperature plasma), 
implies a large value for Tgi- Once Tei becomes comparable or larger than 
the inflow time, electrons and ions can no longer maintain equilibrium 
via Coulomb collisions. Instead, a two-temperature plasma develops in 
which the ions stay close to the virial temperature while the electrons are 
substantially cooler. All current models for low m accretion flows onto 
black holes are based on the two-temperature idea (e.g., Ichimaru 1977; 
Rees et al. 1982; Narayan & Yi 1994; Abramowicz et al. 1995; Bland- 
ford & Begelman 1999; and the numerous variants on these models), and 
would not be viable otherwise. In particular, any collective mechanism 
that coupled the ions to the electrons substantially more efficiently than 
Coulomb collisions would undercut the foundations of two-tcmpcrature 
plasma accretion models. Although no such mechanism is known, the 
well known complexities of laboratory plasmas suggest that the possible 
existence of more efficient coupling cannot be entirely ruled out (Begel- 
man Sz Chiuch 1988; Bisnovatyi-Kogan Sz Lovelace 1997; Quataert 1998; 
Quataert & Gruzinov 1999; Medvedev 2000). 

Although from a purely theoretical perspective the question of whether 
low rh accretion flows form a two-temperature plasma remains open, 
observations — particularly those that indicate a remarkably low luminos- 
ity for the black hole at the Galactic Center (Narayan et al. 1998) — are 
hard to reconcile with almost any imaginable single temperature disk 
model (see below). This has prompted most theorists to accept that 
accretion flows with M <C MEdd develop a two-temperature structure, 
though the exact threshold below which a hot two-temperature flow sup- 



Figure 3.4- Structure of a radiatively inefficient accretion flow, from unpublished ax- 
isymmetric simulations by the author and Kees DuUemond. The numerical method- 
ology is identical to that used by Stone, Pringle, & Begelman (1999). Such flows are 
geometrically thick, turbulent, and weakly bound. Outflows can readily be driven off 
the surface of the disk and along the polar funnel regions. 

plants a cool thin disk is less certain. Narayan & Yi (1995b) suggest that 
the critical accretion rate is around 

"Xcrit ^ ; (3.26) 

that is, one to two orders of magnitude below the Eddington limit. 

Once a two-temperature flow has developed, three model-independent 
conclusions follow. First, since the ion temperature Tj remains compa- 
rable to the virial temperature, the flow can be geometrically thick, with 
h/r ~ 1. Second, since most of the accretion energy remains locked up 
in the ions, the flow is underluminous, with an accretion efficiency much 
smaller than the e ~ 0.1 values realized for thin disks. Finally, since 
the gas is both hot and rapidly rotating, it is at most marginally bound 
to the black hole (formally, the Bernoulli constant, which is conserved 
for adiabatic flows, is positive). As a consequence, it is certainly plausi- 
ble, and possibly inevitable, that a substantial fraction of the putative 
fuel for the black hole actually ends up outflowing from the galactic 
nucleus (Narayan & Yi 1995a; Blandford & Begelman 1999, 2004). Out- 
flows further reduce the already low radiative efficiency (now defined as 
e = L/c^M, where the inflow rate is measured at large radius), because 
most of the gas then fails to sample the deep potential well close to the 
black hole. Figure 3.4 shows semi-quantitatively (not quantitatively, be- 
cause the simulation depicted was two-dimensional and non-magnetic) 
what the resulting flow looks like. 



Apart from such generalities, no consensus has yet been reached as 
to the detailed structure of non-radiative accretion flows (also described 
as advection dominated flows, since most of the liberated gravitational 
potential energy is advected with the fluid rather than being radiated). 
The problem lies in the fact that the structure of non-radiative accre- 
tion flows is expected to depend sensitively on the nature of energy and 
angular momentum transport. The numerical simulations needed to di- 
rectly attack this issue are extraordinarily diflacult. For an illustration 
of the difficulties, suppose that the black hole is accreting from the hot 
interstellar medium of an elliptical galaxy with Tism ~ 10^ K and sound 
speed Cg a few hundred kilometers per second. The ratio of the Bondi 
radius, 

Rb = ^, (3.27) 



to the gravitational radius of the hole — both of which must be resolved if 
we want to get the boundary conditions strictly correct — is then Rs/Rg ^ 
10^, vastly exceeding the dynamical range accessible to any present 
three-dimensional simulation. Two-dimensional simulations can do bet- 
ter but have their own limitations due to both the antidynamo theorem 
in axisymmetry (Cowling 1934) and the fact that the MRI has differ- 
ent properties in axisymmetry, as compared to three dimensions (Good- 
man & Xu 1994; Hawley et al. 1995). Recent examples of simulations 
include work by Stone &: Pringle (2001), Machida, Matsumoto, &l Mi- 
neshige (2001), Hawley & Balbus (2002), Proga & Begelman (2003), Pen, 
Matzncr, & Wong (2003), and Igumenshchev, Narayan, &; Abramowicz 
(2003). There is substantial agreement that outflows of various kinds 
occur generically as a side effect of non-radiative accretion, but opinions 
still differ on the steady-state radial structure of the flow, and on how 
to interpret the numerical results in terms of simpler analytic models. 

Accretion rates above rhcrit (Eq. 3.26) yield flows that are dense 
enough to provide good coupling between ions and electrons. Provided 
that the accretion rate is not extremely high, i.e., m ~ 1, the flow is then 
radiativcly efficient and cools to form a geometrically thin {h/r <C 1) 
disk. The disk extends from an outer radius, which is probably set by 
the considerations of gravitational stability described in the previous 
section, down to the marginally stable circular orbit at r = r^^- Interior 
to the marginally stable circular orbit, which lies at = 6 GMji? 
for a Schwarzschild black hole, circular particle orbits are unstable — any 
inward radial peturbation suffices to place the particle onto a plunging 
trajectory into the black hole. In the simple — but possibly incorrect — 
case where there is zero torque at rms (non-zero torque is discussed in 
§3.4), the radiative efficiency is determined by the binding energy ^^ms 
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Figure 3.5. (Left panel) Location of the marginally stable circular orbit as a function 
of dimensionless spin parameter a for prograde (solid line) and retrograde (dashed 
line) orbits. (Right panel) Efficiency of energy extraction if gas spirals slowly into 
Tms before plunging silently into the black hole. Plotted is 1 — -Ems, whore Ems is 
the binding energy of the marginally stable circular orbit for prograde and retrograde 
motion. 



of gas at that radius. Plots of these quantities are shown as a function 
of the dimensionless spin parameter (a = Jc/GM^ for a hole with an 
angular momentum J) in Figure 3.5. 

The basic theory of thin disks is well known (for reviews, see, e.g., 
Pringlc 1981; Frank et al. 2002). The main open issues (at least for 
flat, circular disks) concern the role of magnetic fields and radiation. 
The structure of magnetic fields near rms is of particular interest, since 
it determines not only the boundary condition at the disk's inner edge 
(Krolik 1999), but also the strength of magnetic fields threading the 
black hole itself. If the black hole is spinning, magnetic fields admit the 
possibility of tapping some of the rotational energy via the Blandford- 
Znajck (1977) process. Spin energy extraction provides one way to power 
jets (e.g., Blandford 2000). and is frequently suggested — mostly on the 
basis of empirical arguments — as the source of the radio- loud/radio-quiet 
dichotomy in AGN (e.g., Wilson Sz Colbert 1995). As a counter to 
such claims, Ghosh & Abramowicz (1997) and Livio, Ogilvie, & Pringle 
(1999) have argued that since dynamo generated disk fields are expected 
to be weak (sub-equipartition) , and dragging in ordered external fields 
is difficult (Lubow et al. 1994), the efficiency of the Blandford-Znajek 
process is likely to be negligibly small for thin disks. These arguments are 
plausible, but I doubt that they represent the last word on the subject, 
not least because of observational indications for spin energy extraction 



in some AGN (Wilms et al. 2001; Reynolds et al. 2004, though note 
that even if spin energy is being extracted, this need not be the classical 
Blandford-Znajek process at work). 

As the accretion rate increases, thin disks are predicted to remain 
the preferred mode of accretion below some upper limit set by the ef- 
fects of radiation pressure. For a spherical flow, that upper limit lies at 
m = 1. In disks, however, matters need not be so simple. Once radiation 
pressure dominates gas pressure in the inner disk, the disk is prone to 
several instabilities, even in the absence of magnetic fields (Shakura k. 
Sunyaev 1976; Agol ct al. 2001). When magnetic fields are included, it 
is found that the flow develops photon bubbles (Arons 1992; Gammie 
1998; Blaes k. Socrates 2001, 2003) — low density radiation filled cavities 
that coexist with dense gas pressure dominated regions. In principle, ra- 
diation can escape from such a highly inhomogenous structure at rates 
above the Eddington limiting luminosity, without being accompanied by 
catastrophic mass loss (Begelman 2001; Begelman 2002; see also Shaviv 
1998). This could permit the existence of stable, thin disks, with an 
emergent luminosity as large as L = 10 — 100 LEdd (Begelman 2002). 
Detailed calculations of the outcome of these instabilities in magnetized 
disks, and their implications for the structure of radiation pressure domi- 
nated flows even below the Eddington limit, are now in progress (Turner, 
Stone, & Sano 2002; Turner et al. 2003; Turner 2004). 

The final regime to consider is that of overfed accretion, with m S> 1. 
By definition, radiation pressure is central to the dynamics of such flows, 
though it need not act to shut off accretion directly. At sufficiently high 
accretion rates, the opacity can be large enough that photons in the 
inner regions are advected inward faster than they can random walk 
away from the black hole. Eventually, they are dragged across the event 
horizon along with the gas and are not radiated from the system. For 
a spherical accretion flow, it is straightforward to show that advection 
wins over diffusion inside a trapping radius given by 

nrap = 7 -77 , (3.28) 

where = 2GM/(? is the Schwarzschild radius (Begelman 1979). The 
same behavior occurs for flows with non-zero angular momentum (Begel- 
man Meier 1982). The emergent luminosity from photon-trapped 
flows remains around L-^dd (e.g., Blondin 1986; Ohsuga et al. 2002), 
so as m increases, the radiative efficiency decreases inversely with the 
accretion rate. 

Clearly, the similarities between the m <C 1 and the rh 3> 1 regimes 
are striking. In both cases, the radiative efficiency is low, for m ^ 1 be- 



cause of inefficient ion-electron coupling, and for m 3> 1 because the high 
opacity prevents photons from escaping the fluid. In the simplest anal- 
ysis, overfed accretion will then also result in a hot, geometrically thick 
disk, with the main difference from the underfed case being the lower 
value of 7 appropriate for a radiation dominated gas. Strong outflows 
could then be inevitable, and these might even reduce the accretion rate 
onto the hole back down to levels comparable to the Eddington limit. 
Once again, simulations of radiation dominated disks will be needed to 
test such speculations. 

In addition to providing an explanation for the dichotomy in the ob- 
served properties of super massive black holes in AGN and quiescent 
galactic nuclei, the existence of these different modes of accretion also 
has implications for the growth of black holes. Accretion rates well below 
the Eddington limiting value probably provide the black hole with a neg- 
ligible amount of mass, despite the fact that this mode lasts much longer 
than more active phases (Haehnelt, Natarajan, &: Rees 1998). More sur- 
prisingly, attempts to flood the black hole with mass at rates vastly in 
excess of the Eddington value may also fail if, as we have speculated, 
such high m flows are as vulnerable to outflows as their low rh counter- 
parts. The entire accretional growth of black holes may be restricted to 
occur via radiatively efficient accretion at rates that lie within an order 
of magnitude or so of the Eddington limit (probably extending above as 
well as below the nominal "limit"). For the most luminous AGN, such 
a "what you see is what you get" model allows for a consistent match 
(within substantial uncertainties) between observations of the luminous 
output of AGN and the local density of supermassive black holes (Barger 
et al. 2001; Merritt k Ferrarese 2001; Yu & Tremaine 2002; following the 
basic approach of Soltan 1982). It may also apply to the earlier growth 
phase of supermassive black holes. 

3.4 Dynamics of the Plunging Region 

The inner boundary of a radiatively efficient accretion disk around a 
black hole lies close to the marginally stable circular orbit. At v ^ ^msi 
the force balance for a fluid element in a thin disk is predominantly be- 
tween gravity and centrifugal force, with angular momentum transport 
driving a slow, subsonic inflow. Interior to rms, circular orbits are un- 
stable. In the standard model for black hole accretion disks (Novikov 
&; Thorne 1973; Page Sz Thorne 1974), once gas reaches this plunging 
region, it flows into the black hole along geodesies, maintaining a con- 
stant energy and angular momentum. In this simple model, gas in the 
plunging region can have observational effects, mainly by reprocessing 
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Figure 3.6. Local flux from the disk (i.e., the emission per unit area in the rest 
frame of the orbiting gas, ignoring returning radiation) for disk models with different 
inner boundary conditions. Lowest curve shows the profile of the emission for a 
standard Novikov-Thorne disk (Page & Thorne 1974) around a Schwarzschild black 
hole, assuming zero torque at rms. F{r) peaks well outside the marginally stable 
circular orbit. The throe upper curves show the effect of adding increasing amounts 
of torque at rms ( Agol & Krolik 2000) , with the labels showing the resulting increase 
Ae in the radiative efficiency as a fraction of the original efHciency e. No emission is 
assumed to originate from the plunging region. 

or reflecting radiation generated elsewhere in the system (Cunningham 
1976; Reynolds & Begelman 1997a). The plunging region does not, how- 
ever, alter the structure of the disk or the radiative efficiency of the flow. 
The inner boundary condition for the disk is simply that there should 
be zero torque at r = rms- 

Early studies of black hole accretion assumed (often explicitly) that 
magnetic fields in the plunging region remained weak. Provided that this 
is so, the inflowing gas rapidly attains supersonic (and super-Alfvenic) 
radial velocities and becomes causally disconnected from the disk out- 
side the marginally stable orbit. This justifies the zero-torque boundary 
condition at rms- If magnetic fields become strong, for example as a con- 
sequence of the shearing of frozen-in fields during the final inspiral, the 
argument fails. Energy and angular momentum can then be exported 




from the plunging region to the disk, modifying its structure and emis- 
sion (Krolik 1999; Gammie 1999a). The same basic process can occur 
in either Schwarzschild or Kerr geometry. 

The consequences of a non-zero torque at the marginaUy stable or- 
bit have been discussed in detail by Agol k. Krolik (2000). The most 
basic is that energy that would otherwise be lost down the hole can be 
transmitted from the plunging region to the disk, and subsequently ra- 
diated. As a result, the radiative efficiency of thin disk accretion could 
be larger than normally assumed (Fig. 3.5). In particular, the standard 
value of e = 0.06 for the yield of accretion onto Schwarzschild holes 
may more usefully be regarded as a lower limit. This implies that — 
although the qualitative trend for more rapidly spinning black holes to 
be more luminous at a given accretion rate than Schwarzschild holes 
is very reasonable — there are significant theoretical uncertainties that 
need to be borne in mind when attempting to determine the average 
spin of black holes from comparison of black hole masses and integrated 
AGN output. For example, the relatively high inferred radiative effi- 
ciency (e > 0.1) derived by Elvis, Risaliti, & Zamorani (2002) and Yu 
&; Lu (2004) need not necessarily imply a rapid average spin for the 
supermassive black hole population. Perhaps more interestingly, given 
that measuring e to a factor of about two precision is obviously a chal- 
lenging task, a non-zero torque at the disk inner edge makes a large 
change to the radial distribution of the disk emissivity. This is plotted 
in Figure 3.6, using expressions given in Agol & Krolik (2000). Unlike a 
"standard" disk, in which the local flux peaks well outside the last stable 
orbit, a torqued disk has a steeply rising emissivity profile right down 
to the marginally stable orbit (at large radius, the contribution to the 
emissivity due to the torque at r^g scales with radius as F{r) oc r~'^/^, 
which is steeper than the dissipation profile for an untorqued disk) . This 
difference in the dissipation with radius is potentially observable. It is 
worth noting that similarly steep dissipation profiles are predicted in 
some models in which magnetic fields act to extract energy from a Kerr 
black hole (Takahashi et al. 1990; Li 2002). 

Numerical simulations have generally been supportive of the basic con- 
cept outlined by Krolik (1999). Pseudo-Newtonian calculations, which 
use a modified potential to mimic some of the most important relativis- 
tic effects (Paczynski & Wiita 1980), have confirmed that there can be 
significant non-zero torques at (Hawley 2000; Reynolds &: Armitage 
2001; Hawley & Krolik 2001, 2002). As shown in Figure 3.7, this leads to 
a declining specific angular momentum profile that extends seamlessly 
through the last stable orbit to smaller radii. The same phenomenon has 
now been seen in General Relativistic MHD simulations of black hole ac- 
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Figure 3.1. (Left •panel) Magnetic torque at the marginally stable orbit as a fraction 
of the torque in the disk at a radius of 2 rms, as derived from numerical simulations 
with relatively high (solid line) and relatively low (dashed line) sound speed (Armitage 
& Reynolds 2003). (Right panel) Specific angular momentum l{r) at several difi^erent 
times for a simulated flow in which the disk magnetic field strength was relatively 
large (Reynolds & Armitage 2001). The inner disk is quite variable, but typically l{r) 
continues to decline interior to rms- 



cretion (De Villiers & Hawley 2003; De Villiers, Hawley, & Krolik 2003). 
Large amplitude (factor ~ 2) fluctuations in the torque also occur and 
persist over a timescale of ~ 10 orbital periods. The importance of this 
torque for the dynamics of the disk may also vary with the tempera- 
ture (or, equivalently, thickness) of the inner disk, with the zero-torque 
limit becoming a better approximation as the thickness decreases (Ar- 
mitage, Reynolds, & Chiang 2001; Reynolds & Armitage 2001; Armitage 
& Reynolds 2003; Afshordi & Paczynski 2003). 

3.5 Variability 

The emission from AGN is variable (Mushotzky, Done, Sz Pounds 
1993; Ulrich, Maraschi, & Urry 1997), as is the much weaker emis- 
sion from the Galactic Center (Baganoff ct al. 2001; Zhao, Bower, & 
Goss 2001; Genzel et al. 2003; Chez et al. 2004). Although in some 
observationally important cases the correlations between variability in 
different wavebands can be understood through simple light travel time 
arguments (Blandford &; McKee 1982), the basic origin of most of this 
variability is not well understood. Here, I discuss what insights into the 
problem recent MHD calculations of disk accretion provide. 

The first question we might ask is, are there large amplitude changes 
in the luminosities of AGN akin to those seen in most other disk accret- 
ing systems? We have already discussed the outbursts of dwarf novae, 



which are attributed to a global thermal instability of the accretion disk. 
Accretion onto some neutron stars and (stellar mass) black holes in low- 
mass X-ray binaries is similarly prone to thermal instability (Tanaka & 
Shibazaki 1996) and can be modeled using a minimally modified version 
of the theory developed for dwarf novae (Dubus, Hameury, & Lasota 
2001). The same ideas have also been pressed into service to explain the 
outbursts of protostellar disks seen in FU Orionis objects (Hartmann & 
Kenyon 1996; Bell &: Lin 1994), though the identification of these events 
with thermal disk instabilities is less compelling. 

If thermal instabilities act on a global scale in AGN, the timescalcs for 
the cycling between outburst and quiescence would be long — probably 
~ 10^ yr or longer (Lin & Shields 1986). Transitions between "on" 
and "off" states might not occur on easily observable timescalcs, though 
they could leave traces in the form of intermittent activity from radio- 
loud AGN (Reynolds & Begelman 1997b; Owsianik, Conway, & Polatidis 
1998). If such intermittent accretion is commonplace, there ought to 
exist a population of relatively quiescent galaxies in which the black 
hole is surrounded not by a radiatively inefficient fiow, but rather by 
a temporarily inactive thin disk. Detailed models of this type have 
been developed by Sicmiginowska, Czerny, &; Kostyunin (1996) and by 
Siemiginowska &: Elvis (1997). 

Disks in AGN extend out to radii that encompass the zone in which hy- 
drogen is partially ionized, so the basic ingredient necessary for thermal 
instabilities to occur is present (Burderi, King, & Szuszkicwicz 1998). 
What is less clear is whether the efficiency of angular momentum trans- 
port (i.e., the value of a) varies between the ionized and neutral state in 
the AGN environment. A substantial difference between ahot and cccoid is 
needed if the whole disk is to flip between outburst and quiescent states. 
In dwarf novae, the required change probably occurs because the hot 
and cold states lie on opposite sides of the critical degree of ionization 
required for sustained MHD turbulence (Eq. 3.22; Gammie & Menou 
1998). It has been argued (Menou & Quataert 2001) that this is no 
longer the case in AGN disks, leading to thermal instabilities that pro- 
duce local flickering but not global outbursts. If so, any observational 
evidence for large scale outbursts in AGN may require a different ex- 
planation. Numerous possibilities for time dependent behavior suggest 
themselves in the outer disk, where there is an interplay of gravitational 
and MHD disk instabilities, though none have been worked out in any 
great detail. 

At frequencies between ~ 10~^ Hz and ~ 10^*^ Hz (i.e., timescalcs 
of tens of minutes to around a year). X-ray observations provide good 
constraints on the properties of AGN variability (e.g., Edelson & Nan- 



Figure 3.8. View of a simulated disk as seen by a distant observer at an inclination 
angle of 55° (Armitage & Reynolds 2003). The patchy emission is boosted on the 
approaching (left hand) side of the disk as a consequence of beaming. In constructing 
this image, the local disk emission has been assumed to scale with the vertically 
integrated magnetic stress. Raytracing through the Schwarzschild metric is then 
used to account for the relevant relativistic effects, including beaming, gravitational 
redshift, and light bending. 

dra 1999; Uttley, McHardy, & Papadakis 2002; Vaughan k Fabian 2003; 
Markowitz et al. 2003; and references therein). The power spectral den- 
sity G{f) derived from these observations is often consistent with a power 
law broken at a frequency /break, 

G{f) OC /< /break, 

G(/) OC /-""-"^ /> /break. (3.29) 

Representative values for the two slopes are aiow — 1 and Ohigh > 2 
(Vaughan & Fabian 2003; Markowitz et al. 2003). 

Although a detailed comparison is beyond the capabilities of current 
numerical simulations, there appear to be good prospects for relating 
noise power spectra of this form to theoretical models of MHD disk 
turbulence. Simulations of magnetically active disks show that the flow 
displays variability across a wide range of timescales, with physical quan- 
tities such as the magnetic stress and mass accretion rate showing tem- 
poral power spectra that are described by power laws (Kawaguchi et al. 



2000; Hawley & Krolik 2002; Armitage k Reynolds 2003). The patchy 
and rapidly fluctuating pattern of stress across the surface of the in- 
ner disk, which gives rise to this variability, is shown in Figure 3.8. If 
we assume that the local disk emission traces the magnetic stress, then 
the power spectra predicted from the simulations have the form shown 
in Figure 3.9. Regardless of the inclination of the disk to the line of 
sight (which affects the importance of relativistic effects such as beam- 
ing), the temporal power spectrum at frequencies comparable to those of 
the inner disk is around G{f) oc f~^. Individual annuli in the simulated 
disk, moreover, produce power spectra that break — at around the orbital 
frequency — to substantially steeper slopes of around /~^'^ (Armitage & 
Reynolds 2003) . These results support the common assumption that the 
break frequency in observed systems scales with the orbital frequency at 
the marginally stable orbit, and hence with the mass of the black hole 
(e.g., Nowak & Chiang 2000; Lee et al. 2000). Although there are signif- 
icant uncertainties in many AGN black hole masses, the observations of 
Markowitz ct al. (2003) suggest a correlation between black hole mass 
and break frequency. We note, however, that work to date has not been 
able to recover the observed values of aiow and awgh — the theoretical 
power spectra are steeper than observed. 

3.6 Warped, Twisted, and Eccentric Disks 

Flat disks, in which the gas follows close to circular orbits, are special 
cases of the more general situation in which the disk may be eccentric 
and/or warped out of a single plane. A common assumption is that an 
initially warped or eccentric disk will relax to a flat circular one "on a 
viscous timescale" . This is at best a gross simplification of the evolution 
of warped or eccentric disks, and may even be qualitatively wrong. As 
one might expect, the evolution of warped or eccentric disks is rather 
complex, and in this section I attempt only to provide a framework for 
understanding some of the key results. For more details, the interested 
reader is well advised to consult the recent and exhaustive studies by 
Ogilvie (1999, 2000, 2001), which include references to the earlier liter- 
ature on the subject. 

Observations provide clear motivation for considering the evolution 
of warped disks in AGN. The masing disks in NGC4258 (Miyoshi et al. 
1995), NGC1068 (Greenhill & Gwinn 1997), and the Circinus galaxy 
(Greenhill et al. 2003) all appear to be warped, albeit to different de- 
grees. This raises two questions. First, are the warps merely decaying 
features that reflect the complex angular momentum distribution of the 
infalling gas that formed the disk, or are they self-excited by some pro- 
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Figure 3.9. Temporal power spectra G(/) derived from a numerical simulation of 
accretion onto Schwarzschild black holes (Armitage & Reynolds 2003). The frequency 
/ has been rescaled such that the orbital frequency at the radius of marginal stability 
corresponds to / = 1. Shown here are the power spectra for i = 1°, i = 30° (offset 
by 10^), i = 50° (offset by 10^), and i = 80° (offset by 10*). MHD turbulence in 
the disk generates variability across a wide range of timescalcs, with quantities such 
as the predicted emission or mass accretion rate showing steep, approximately power 
law spectra. 

cess intrinsic to the disk or AGN itself? Second, if the warps are in fact 
decaying, how rapidly does that process occur compared to the timescale 
on which the gas would be accreted? 

Although easily formulated, both of these questions have proved to 
be difficult to answer. The basic physics revolves around the fact that in 
a warped disk, neighboring annuli are inclined relative to one another, 
and as a consequence, there is shear in the vertical direction as well as in 
the radial direction. Dissipation of the energy associated with the out- 
of-plane motions is governed by the (r, z) component of the stress, and 
this acts to flatten out the disk as, simultaneously, the (r, 4>) stress drives 
accretion. The crucial point is that these stresses can be generated by 
entirely different physical processes (for a concrete example, see Gammie, 
Goodman, & Ogilvie 2000) and need not act on the same timescale, even 
roughly (Pringle 1992). The detailed physics of stresses within the disk 
will then determine whether the disk flattens before a significant fraction 




of the gas accretes, or whether instead the warp is advected inward along 
with the gas. 

To go further, we need to elaborate on the mechanisms that can al- 
low annuli to exchange angular momentum. If the disk is thin enough, 
specifically, if 

-<a, (3.30) 
r 

then the evolution of the tilt as a function of radius, like the evolu- 
tion of the surface density, is diffusive (Papaloizou & Pringle 1983). If, 
conversely, 

->a, (3.31) 
r 

then the warp is communicated through the disk by waves (Papaloizou &: 
Lin 1995). Assuming that a ~ 0.1 in geometrically thin AGN disks, then 
most parts of the disk (except, perhaps, the innermost regions) are likely 
to fall into the diffusive regime, which we will consider exclusively from 
now on. We will also ignore the complications of self-gravity (Papaloizou, 
Terquem, &; Lin 1998), even though, as already noted, the masing disks 
cited as motivation are probably self-gravitating. 

In the diffusive limit, the evolution of the disk can be described to a 
first approximation using the simple theory described by Pringle (1992). 
In this theory, the stresses within the disk are reduced to two viscosities; 
ui, which describes the usual kinematic viscosity leading to inflow, and 
f2, which is related to the stresses which act to flatten out any warp. The 
timescales for accretion t^, and for decay of the warp twarp are related to 
these viscosities in the usual way: 

U = (3.32) 

.2 



^warp 



r 

1^2 



The relation between these two viscosities can be calculated either an- 
alytically (Ogilvic 1999) or via numerical simulations (Torkelsson et al. 
2000). Both approaches suggest that 

^=.^»l. (3.33) 



This is a striking result. It implies that warps are likely to decay rapidly 
and will require strong forcing if they are to persist for long periods. 
Torques from the reprocessing of radiation from the central source by 



the disk provide one well studied forcing term, which could act to ex- 
cite warps even in initially planar disks (Pringle 1996; Maloncy, Bcgel- 
man, & Pringle 1996). The large ratio of 1^2 /'^i implied by the above 
equation, however, severely limits the circumstances in which radiation 
driven warping can overcome the disk's intrinsic tendency to flatten out. 
Torques from a disk wind are potentially more promising (Schandl &; 
Meyer 1994; Lai 2003), though harder to evaluate in detail. 

Further complications ensue if the black hole is rotating, even mod- 
estly. Orbits whose angular momentum is inclined to the spin axis of a 
rotating hole will process due to the Lense-Thirring effect ( "frame drag- 
ging"). This forced precession is only significant at small radii, where 
it will affect the inclination of the disk. Again, the simplest case to 
analyze is the one in which the inner disk is in the diffusive regime of 
warp evolution. Differential precession, acting on an initially flat but 
misaligned disk, will rapidly induce a sharp twist in the disk close to the 
marginally stable orbit. This will be flattened out by viscosity, with the 
consequence that the disk near a spinning hole will be driven into the 
equatorial plane of the hole, regardless of the inclination at large radii 
(Bardeen k, Petterson 1975). Over longer timescales, the torque between 
the disk and the hole will also act to realign the system until the angular 
momentum vector of the hole, and of the disk at all radii, are coinci- 
dent. Since normally Jdisk 3> J,, this balancing act will usually result in 
a large change in the spin axis of the black hole until it roughly matches 
the angular momentum vector of gas in the outer disk. How quickly this 
process occurs depends, again, largely on ^2/1^1, with large values of this 
ratio leading to rapid alignment (Natarajan & Pringle 1998). 

In some circumstances, the inner accretion disk around a rotating 
black hole may be thick enough that warps display wave-like evolution. 
In this case, the inner disk may be able to remain misaligned with respect 
to the spin axis of the black hole, even in the presence of dissipation 
(Demianski & Ivanov 1997; Ivanov & Illarionov 1997; Lubow, Ogilvie, k. 
Pringle 2002) . This could have important consequences for the direction 
of jets launched from the inner regions of the disk. 

Unlike in the case of warped disks, there is little observational ev- 
idence to suggest that AGN disks — or any disk that is not strongly 
forced, for example by a binary companion — are significantly eccentric, 
though some sort of asymmetry is deduced in as many as 60% of AGN 
with double-peaked Balmer lines (Strateva et al. 2003). Theoretical 
studies, however, raise the interesting possibility that there could be 
circumstances in which isolated accretion disks spontaneously develop 
eccentricity. In particular, a two-dimensional disk model in which the 
viscosity follows a Navier-Stokes form, with shear viscosity jx and bulk 



viscosity is unstable to the growth of eccentricity if 
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that is, if the vertically integrated viscosity is increasing too rapidly with 
surface density (Ogilvie 2001). If the bulk viscosity is small, as is often 
assumed in accretion disk modeling, then this result implies that almost 
all disks ought to be unstable to the growth of eccentricity (Lyubarskij, 
Postnov, k Prokhorov 1994; Ogilvie 2001). 

Earlier, of course, we specifically warned of the dangers of assuming 
that angular momentum transport in an MHD flow can be described 
via an effective Navier-Stokes viscosity. Indeed, Ogilvie (2001) shows 
that if angular momentum transport follows an alternative analytic pre- 
scription, which is motivated by the phenomenology of MHD turbulence, 
then the wholesale instability of disks to eccentricity can be avoided. The 
important point is that the behavior of eccentric disks — which at a min- 
imum will exist in galactic nuclei following the tidal disruption of stars 
(Rees 1988), or when binary black holes arc present — is closely related 
to subtle aspects of turbulent transport within the disk. Predictions as 
to their evolution cannot yet be made with confidence. 
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